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FOREWORD 


This computer program description is submitted by the AiResearch 
Manufacturing Company, Los Angeles, California, a Division of the Garrett 
Corporation in compliance with Part IV; paragraph A of NASA Contract NASl- 
5116 and paragraphs 6.4.6 and 5.4.1 of NASA Statement of Work L-4947. 

This report is one in a series of several computer program description 
reports that will be forwarded to NASA in compliance with contractual 
requirements. This report will be identified at AiResearch by AiResearch 
Report Number Ap-1001-t. 
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hypersonic ramjet experiment project 

PHASE I 

COMPUTER PROGRAM DESCRIPTION 
RAMJET AND SCRAMJET CYCLE PERFORMANCE 

By: Robert J. Jackson 

Tennyson T. Wong 

SUMMARY 


The Ramjet and Scramjet Cycle Performance Program performs one- 
dimensional calculations of the equilibrium, real-gas internal flow properties 
and of pertinent engine performance parameters. This program is a modified 
version of the NASA program developed by L. Franc iscus. 

The program can be used (1) for preliminary design calculations, i.e., 
the sizing of critical flow areas in the engine to comply with stipulated 
internal flow properties and (2) for design analysis, i.e., determination 
of internal flow properties corresponding to stipulated flow areas. Only the 
combustion of hydrogen in air is considered. 
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PROBLEM DESCRIPTION 


General capabilities and input-data requirements of the program are dis 
cussed in the following paragraphs in terms of the major engine components: 
(I) inlet, (2) combustor, and (3) exhaust nozzle. Pertinent engine stations 
are shown below. 





Two-Dimensional Schematic of Engine Stations 

Inlet calculations .- Flight conditions are specified in terms of flight 
Mach number, ambient pressure, and ambient temperature. 

As shown above, station 2 is the interface between the engine inlet and 
the combustor. The program calculates the one-dimensional, equilibrium prop- 
erties of the inlet air at station 2 corresponding to 0) assigned flow area 
or Mach number, (2) inlet losses, and (3) inlet spillage (mass-flow ratio). 
Note that station 2 is primarily the combustor inlet, not necessarily the 
engine-inlet throat (minimum flow area). 

The inlet losses are stipulated in terms of the kinetic energy effi- 
ciency, 1W, or the rise in enthalpy, AH’, from conditions corresponding to 
free-stream pressure and entropy to conditions corresponding to free-stream 
pressure and combus tor- i n let entropy. 

The inlet spillage is given as the ratio of capture area to cowl area. 

Data printed out at station 2 include pressure, temperature, enthalpy, 
entropy, molecular weight, Cp, y, Mach number, ratio of flow area to inlet 
capture area, velocity, total pressure recovery, vacuum specif ic impulse, anc 
fuel port area. Also, the equilibrium Mol fractions of the working-fluid 
constituents are printed. 


2 



Combustor calculations . - The combustor calculations yield the one* 
dimensional, equilibrium properties of the combus tor^ex i t flow, engine 
.station 3. The following conditions and data must be input: 

a. Mode of combustion, supersonic or subsonic. 

b. Fuel equivalence ratio. 

c« Fuel storage temperature and fuel temperature at injection ports. 

d. Fuel injection velocity. 

e. Ratio of combustor wall area to inlet flow area. 

f. Combustor wa I l-f r i ot ion drag coefficient. 

g. Assigned combustor exit condition (pressure, Mach number, or area). 

h. Combustion chemical efficiency, Ihe ratio of reacted fuel to fuel 
supp 1 ied. 

The input provision for two values of fuel temperature is necessary for 
analysis of regenerative ly-cnoled combustors. For energy calculations, the 
fuel enthalpy must be based on the temperature of the fuel prior to its 
introduction into the cooling circuits; however, sizing of fuel injection 
ports must be based on the actual fuel injection temperature and velocity. 

The equilibrium thermodynamic properties of the combustion products are 
calculated by the procedure described in ref. I. 

Calculation of combustor exit momentum requires, for cases other than 
constant-area or constant-pressure combustion, evaluation of combustor-wall 
pressure forces. In general, this is done by assuming that the effective 
combustor static pressure is the average of the combustor inlet and discharge 
static pressures. However, where appropriate, the effective combustor static 
pressure can be taken to be: 

p _ P 2 + KP 3 • (I) 

s ' K + I 


By formula ( I ), the effective combustor static pressure can be made to 
approach the combustor exit pressure, Pj, as K is allowed to increase. Use 
of this procedure would be based on detailed investigation of the combustor 
pressure distribution. 


The combustor wall-friction drag coefficient, Cq, is defined by: 



A 


w 


q 


( 2 ) 
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where F n is the wall-friction drag force, A* the wall area, and the average 
of combustor inlet and discharge dynamic pressures. 


The flow properties printed out at 
printed for station 2 , the area ratio, 
eters: Net internal thrust, thrust per 

impulse, and thrust coefficient. 


station 3 include the parameters 
A 3 /A 2 , plus engine performance pa ram 
unit of capture area, fuel specific 


N nzzle calculations .- An isentropic, equilibrium, real-gas expansion 
from the combustor exit to the nozzle exit is calculated. For subsonic com- 
bustion, the first line of output subsequent to the ccmbustor exit refers to 
the nozzle sonic throat (this is omitted for supersonic combustion). Subse 
quent lines of output correspond to assigned values of combustor ex t 
nozzle pressure ratio, to full expansion to ambient pressure, and final ly, 
the assigned nozzle exit condition (assigned area or pressure). For super- 
sonic combustion, the ratio of nozzle area to combustor exit area P r ' n • 
for subsonic combustion, the ratio of nozzle area to nozzle sonic throat 

printed out. 

Nozzle losses are accounted for by a velocity or stream-thrust coeffi- 
cient, C^: 


c s * 


V , RT 

Kac = a V - 

Iva c,ID (^ + ^)id 


( 3 ) 


where the ideal values are calculated for isentropic, equilibrium nozzle 
expansion. The coefficient, C s , is evaluated from detailed nozzle investi- 
gations to account for kinetic, friction, and divergence losses. 

To calculate the net thrust for flight conditions where inlet spillage 
must be considered, an additive drag coefficient, C DA , is employed. This 
drag coefficient is defined as follows. 


F DA ~ C DA q o A c 


( 4 ) 


where Fp A Is 


the drag force caused by inlet spillage, q Q is the free-stream 

dynamic pressure, and Ap is the cowl area. Values of Cq A have been 
correlated with flight Mach number and inlet mass-flow ratio. 


OPERATING INSTRUCTIONS 


Hardware requirements . - 

a. IBM 7090 or 7094 computer with 32K core. 

b. Data Channels and tape units for the software. 
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No on-line printer and punch. 


c » 


d. Peripheral equipment for card-to-tape* tape-to-pr i nter * and tape-to 
punch as required by standard IBM system* 

e. It is assumed that standard keypunch and verifier machines* card 
readers* and printers are used; therefore* these have not been 
specified. 


Software reau i remen ts . - 

a. IBSYS Operating System Tape Version 12 

(1) Basic Monitor ( IBSYS ) Version 4 

(2) Processor (IBJ0B) Version 4 

(3) Assembly Language (MAP) Version 4 

(4) Assembly language (FTC) Version 4 

(5) Loader (ifiiDR) Version 4 

b. The Fortran IV I/O Library Subroutines 

Tape usage .- The standard Fortran input/output configuration* as 
described in IBM 7Q90/7Q94 IBSYS Operating System : IBJOB Processor . Form 

C28-6275* is used for the logical tape unit designation and the file specifi- 
cations are as follows. 


UNITOI 

FILE, 

UTI. 

READY, 

INOUT, 

BLK 


2 56, 

BIN 


UNIT02 

FILE, 

UT?, 

READY, 

INOUT , 

BLK 

= 

256, 

BIN 


UNIT03 

FILE, 

UT3, 

READY, 

INOUT, 

BLK 

= 

256, 

BIN, NOLIST 


UNIT04 

FILE, 

UT4 , 

READY, 

INOUT, 

BLK 


2 56, 

BIN, NOLIST 


UNIT05 

FILE, 

IN, 

READY, 

INPUT, 

BLK 


14 , 

MULTIREEL, BCD, 

NOLIST 

UNIT06 

FILE, 

OU, 

READY, 

OUTPUT, 

BLK 

= 

1 10, 

MULTIREEL, BCD, 

NOLIST 

UNIT07 

FILE, 

PP, 

READY, 

OUTPUT, 

BLK 

= 

28, 

MULTIREEL, BIN, 

NO LIST 

UN IT 08 

FILE, 

* 

MOUNT, 

INOUT, 

BLK 

— 

22, 

BCD 



The logical tape unit 5 and 6 are used as input and output* respectively. 
No intermediate storage tape is used . 
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Deck set-up.- 


a. $J£)B 


b. $EXECUTE 

IBJ0B 

$ IBJ0B 

t 

(Opt ions ) 

1 

Program 


decks 

I 


▼ 

End of File 

* 

(7,8 Punch in Col . 1 ) 

1 

Data 

/Thermodynamic data and\ 

decks 

1 

Vprogram input data / 

i 

End of File 

(Same as above) 


For a description of these control cards see: 

(IBM 7090/7094 Operation Systems Basic Monitor (IBSYS), 
Form C28-6248) 


c. 


(IBM 7090/7094 IBSYS Operating System IBJtfB Processor, 
Form C28-6275) 


$IBSYS 
d. $ST0P 



Deck Set-up 


* PROGRAM DECK CONSISTS 
OF MAIN PROGRAM MAIN 
AND FIVE SUBROUTINES. 

CflREZ 

MATRIX 

GAUSS 

aTeF * 

COMPOS 
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Timing .- The computer time for each case depends upon the number of 
iterations. It takes approximately one minute per case. 


Multiple cases .- The cases may be stacked for one computer run. 


Input Descr ipt ion 


Format . - 
Card I 
I 

914 


Field 


1 

2 

3 

4 

5 

6 


Descr ipt ion 

Reference case number (usually "I”) 

"0" means subsonic combustion; " I means 
supersonic combustion. 

B lank 

Blank 

Number of nozzle stations for which data are 
printed out. 

Number (0 to 25) of nozzle pressure ratios for 
schedule of assigned nozzle pressure ratios. 


7 Blank 

8 "0" means normal printout. "I" means 

"debug printout". 

9 Number of cases (l to 15.). 


2 

8FI0.0 


I 


2 

3 


"0." means specified combustor pressure ratio. 
"-I." means specified combustor exit Mach 
number. "I.” means specified combustor area 
rat io. 


Ae/AC, ratio of nozzle exit area to cowl area. 
Blank if nozzle pressure ratio is specified. 


A2/AC, ratio of combustor inlet area to cowl 
area. Blank if combustor inlet Mach number is 
specif ied . 
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Card Field Description 


2 

8FI0.0 

4 

5 

M2, combustor inlet Mach number. Blank if 
A2/AC is specified. 

Pe/PO, ratio of nozzle exit pressure to 
ambient pressure. Blank if A /AC is 
specif ied. 

3 

1 

cos of fuel injection angle 

8FIO.O 

2 

AC, cowl area, sq ft 

4 

1 

Aw/A2, ratio of combustor wall area to com- 
bustor inlet area. 

5 

1 

Ambient pressure, atm. 

8FI0.0 

2 

Ambient temperature, °K 


3 

Flight Mach number 


4 

B lank 


5 

Additive drag coefficient 

6 

1 

TIKE, inlet kinetic energy efficiency, or 
(by pre-fixing a minus sign) the enthalpy 

8FI0.0 


rise between free-stream pressure, entropy 
conditions and conditions corresponding to 
free-stream pressure and combustor- in let 
entropy. 


2 

Ao/AC, ratio of inlet capture area to cowl 
area ( i . e . , mass-flow ratio). 


3 

Ae, /AC, f , ratio of geometric nozzle 

exi t? e a?ea to effective nozzle exit area (to 
account for boundary- layer growth). 


4 

C s , nozzle velocity coeff i cient, rat io of I yac 
to I...,, calculated with no nozzle losses. 

VdC 
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Card 


Field 


Descr i pt i on 

6 5 Temperature of fuel at fuel -inject ion 

nozzles, °K. 


7 

8FI0.0 


Cards 8-10 
(Optional ) 
8FI0.0 


1 0, equivalence ratio 

2 T) c , combustor chemical efficiency, ratio (by 

weight) of reacted fuel to fuel supplied. 

3 Vf, fuel inlet velocity, ft per sec. 

4 Temperature of fuel corresponding to storage 

condi t ions. 

5 Assigned combustor pressure ratio, area ratio, 
or exit Mach number, as applicable. 

6 Combustor wal 1 -fr i ct i on drag coefficient. 

Up to 25 values of nozzle pressure ratio, 

Used when sixth field of card I is non-zero. 


The "reference case" comprises a complete set of the above seven data 
cards. Up to 14 additional cases, representing changes in any or all of the 
variables entered on Cards 5, 6, and 7, can be handled for constant values of 
the data on cards I to 4. For each of these additional cases, a set of cards 
5, 6, and 7 is added to the input-card deck. However, only the values which 
differ from those of the reference case are punched. For example, to run a 
range of M3, two blank cards and a card punched with the appropriate value in 
the fifth field are added to the input-card deck for each new value of M3. 

Any number of the above i nput-data-card decks can be stacked for a com- 
puter run. 

Table I presents a computer input data format with program internal 
symbols. 

Def ini tions . - 

Additive Drag = C qA (l/2 p^ 2 ) A c 

Combustor Wall-Friction Drag = C , q A 

3 a w 

where q = 1/2 [(1/2 * (1/2 P^ 2 )] 

A = Combustor internal wetted area 
w 
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TABLE I 

IBM COMPUTER PROGRAM INPUT FORM 


IBM MACHINE rm,c r 

DATE. P « 0GR ' 

PREPARED BY PERM- N0 - EXT - . 

EST. RUN TIME CHARGE NO. 

COMMENTS INPUT DATA F0RHAT WITH PROGRAM INTERNAL SYMBOLS 


PAGE NO. _ 
PROGRAM NO. 
EXT. 


1 to 

I < 20 

21 30 

31 40 

41 50 

51 60 C 

! 70 j 

71 80 



*> + •J’ 

s* / 

/ 

4? 

Jr 

❖ 

J 

\ F0RMA1 

r 9t4 

C0NST 

ZEIAR 

A2AC 

VMACH2 

hpam 




C0STH 

AC0WL 







AW? 







- 

(PCC/, 

(TCC ) i 

(VMACH0), 


(CD0)i 




(ANKE), 

(A0AC )t 

(AEA*), 

(CVN)i 

(TTFUEL), 




(PHI), 

(DELTAH), 

(VF), 

(TFUEl), 

(P3P2). 

(CD), 



. 

. 


- 



► FORMA' 

r 8FI0.0 


* 

• 

- 





* 

. 

. 

* 

* 




(pccv* 

(tcc) n 

(vmachU^ 


(COB^ 




(anke) n 

(abac ) n 

(AEA4) m 

(cvn) h 

(TTFUEL^ 




(phi) n 

(deltah) n 

<">» 

(tfuel) n 

(P3P2) n 

(cd) m 


— 


WHERE Nr ICASES 


K..T) AIRES! ARCH MftNI II At IIjRlNl. DIVISION 
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Output Description 


Input data .- The input data values for a run are printed out prior to 
and separate from the calculated results. 

Inlet calculations .- First line (below headings): 

PO - Free stream pressure, atm 
TO - Free stream temperature, °K 
ETA KE - T1 

MACHO - Flight Mach number 

MACH2 - Assigned combustor inlet Mach number (if applicable) 

PHI - Equivalence ratio 

A2/AC - Assigned ratio of combustor inlet area to cowl area (if 
appl i cable) 

VO - FI ight veloci ty 

HOI - Enthalpy (cal per gm) corresponding to free stream pres- 
sure and combustor inlet entropy. 

AO/AC - Ratio of capture area to cowl area. 

AO - Capture area, sq ft 

The second and subsequent lines of printout (excluding the composition 
data) are the values calculated for each iteration to the applicable com- 
bustor inlet condition. The converged values are presented in the last such 
line of data. Units of enthalpy and entropy are cal per gm and cal ger gm -°K, 
respectively, where enthalpy is assigned the value of zero at 298.15 K. 

PO/P - Ratio of free-stream pressure to combustor inlet pressure. 

P ATM - Combustor inlet pressure, atm 

MACHNO - Combustor inlet Mach number 

A/AO - Ratio of combustor inlet area to capture area 
The composition data are the mol fractions of each constituent, 

1 1 



Finally, the inlet total pressure recovery is the calculated ratio of 
combustor inlet total pressure to the free-stream total pressure; fuel inlet 
port area is given in sq ft; and IVAC is the momentum parameter, vacuum 


specific 


impulse. 


V PA 
g W' 


sec. 


Combustor exit calculations .- The heading stipulates the type of nozzle 
calculation; the second line gives the amount of fuel in terms of three 
parameters: (I) the air-to-fuel mass ratio, (2) the equivalence ratio, and 
(3) the f ue 1 -to-m?xture mass ratio. 


Combustor exit data are printed in two lines. Parameters not already 
identified are the following: 

A/A2 - Ratio of combustor exit area to combustor inlet area. 
T/AO - Net thrust per unit of capture area, lb per sq ft. 

SP IMP - Net thrust per unit of fuel flow, sec 
THRUST - Net thrust, lb 


IVAC 

CT 


Vacuum specific impulse, sec 

THRUST 


- Thrust coefficient. 


q A 

■/"i 


i '“'r 

O C 


The second line of output (following composition data), in the case of 
subsonic combustion, is the nozzle throat calculation. Subsequent lines of 
output comprise (I) stipulated nozzle pressure ratios (per word 5 on card I) 
and (2) results of each iteration to the specified nozzle exit condition. 

For subsonic combustion, the area ratios, A/A3, are actually ratios of nozzle 
area to throat area. In the case of supersonic combustion, there is no nozzle 
throat calculation, and A/A3 is actually the ratio of nozzle area to com- 
bustor exit area. 


Test Case 


The following test cases are used to illustrate the input and the output 
of the program: 

a. Subsonic combustion with assigned combustor exit Mach number. 

b. Supersonic combustion with assigned combustor exit Mach number. 

c. Supersonic combustion with assigned combustor area ratio. 



The IBM input forms for these cases are presented in tables II, HI, and 
IV. 


The outputs of the test cases are forwarded under separate cover 

•v 


REFERENCE 


I. Zeleznik, Frank J., and Fordon, Sanford, "A General IBM 704 or 7090 
Computer Program for Computation of Chemical Equilibrium Composition, 
Rocket Performance, and Chapman-Jouge t Detonations", NASA TN D-1454 
( 1 962 ) . 


December, 1965 

Ai Research Manufactur i ng Company 
Division of the Garrett Corporation 
Los Angeles, California 
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TABLE II 

IBM COMPUTER PROGRAM INPUT FORM 
TEST CASE NO. I 


IBM MACHINE PAGE NO. 

DATE PROGRAM NO. _____ 

PREPARED BY PERM. NO. EXT. 

EST. RUN TIME CHARGE NO. 

COMMENTS SUBSONIC COMBUSTION WITH ASSIGNED COMBUSTOR EXIT MACH NUMBER 





























TABLE III 

IBM COMPUTER PROGRAM INPUT FORM 
TEST CASE NO. 2 


IBM MACHINE 


PAGE NO. 

DATE 


PROGRAM NO. 

PREPARED BY 

PERM. NO. 

EXT. 

EST. RUN TIME 

CHARGE NO. 



COMMENTS SUPERSONIC COMBUSTION WITH ASSIGNED COMBUSTION EXIT MACH NUMBER 


1 10 

1 1 20 

2 1 30 

31 40 

'■ ! 50 

5 1 60 

61 70 

71 80 

1 1 

1 0 

4 0 

0 0 

2 


1 

J 

► FORMAT 914 

-1 . 

2.0 

.IU 







0 . 

1,76 








35.0 








’ 

.0334 

220. 

6.0 

0 . 

, 025 



> FORMAT 8FI0.0 

.96 

.82 

1.0 

,9696 

' 

556. 





.6 

.97 

54 50. 

20. 

1 .2 

.002 






BLANK 

CARD 








BLANK 

CARD 



J 



1 , 



■■■■ r 

















































t 

i 
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TABLE IV 




IBM COMPUTER PROGRAM INPUT 

FORM 



TEST CASE NO. 3 


IBM MACHINE 


PAGE NO. 

DATE 



PROGRAM NO. _ 

PREPARED 

BY 

PERM. NO. 

EXT. 

EST. RUN 

TIME 

CHARGE NO. 


COMMENTS 

SUPERSONIC 

COMBUSTION WITH ASSIGNED COMBUSTOR AREA RATIO 


.955 I .0 


1.0 


20 ?! 30 


0 


4 0 I A I 50 


3 


0. 


.93 


2000.0 20.0 


blank! card 


BLANK CARD 


blank[ CARD 
CARD 


FORMAT 914 


FORMAT 8FI0.0 


APRf SEARCH MANUFACTURING HI VISION 

1 .V, rjUrv"* 



































Reference: (l) Aiftssearch Report AP- ICO ) - 1 , Computer 

pi*Op roift Description, Ramjet end Scrcmjot 
Cycle Performance conducted under NASA 
Contract KAS 1 -5 1 1«$ 

Revised Cycle Performance Program (from Reforonco one) including: 


S o) Computer tape 

b) program listing 

c) Sampi o calculation 

d) Input format/ instruct ions 
e$ Sample Input cords 


The revised cycle performance from the referenced report Incorporates 
the following mod if leations; 

l. - The program was modified in order to more ««urately 

for the pressure distribution along the mu It. -staged colter 
walls This can nw be accomplished by stnc,;»ng the inpuv 

doll hr cid. .us., SKCiftino «* ***»>«'. sres r«,. «.4 

ecu ivc lent pressure, and weighing constant * ^ 

Tha orocram will execute the muW i-stajo prc>al*m a •>. 
1514 by «1.«0 th. results fro,, previous stoics. 


2, The momentum-continuity equation in this program was 
modified to a third order algebraic equation in order 
to moro accurately solve for tho velocity at the cost- 

. bustor exit. 

3, Because tho dynamic pressure distribution Is not linear 
along the multi-stage combustor, the viscous drag cal- 
culation was modified by using the averaged vaiuo of q 
o_f_ coph-S.tago rather than using the averaged value of 

q for the overall combustor os dona previously. 



p {"M 3 ? * : •* \ \r 

t. i $ i - 1 - ‘ + . i ^ t 


! , •■■" ; i 

1 jS; ! 

«-/ « V * J • j Z?) **♦ ?»L 
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| vrv 


AIR 


o 

J 


■A 1 '“] r t1 r <- r r r ■ , , „ , 
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C5 A M li E T T C D n p fJ K ATiON 


(* Ci i - & 0 i> i stPu l v Ti n a ii a u l r V/» <-» d • i*m. a ro c- q l r r; . calipohnia n n a u o 
TCLIPHDN t: I21C#| 7?0*10 lO. 670*0131 ■ C t iU F : S A fl M [ T T A I P LOB ANGELES 


In reply refer to* 
SftJD~8267~0510 

10 Hoy 1968 


Motional Aeronautics and Space Administration 
Langley Research Center 
Hampton, Virginia 23365 


Attention: Hr. it. 0. Albert 

Con 1: re c t Adm 1 n i s t ra tor 

Subject: Hypersonic Research Engine Project 

NASA Contract MAS 1-6666 
Submittal of Computer Tope 

Reference: (l) A 1 Research Report AP~lC0l-l, Computer 

P ro grain Description, Ramjet and Sc ramjet 
Cycle Performance conducted under MAS A 
Contract HAS l -5 11 6 

Gentlemen: 

Tho following Information is submitted as requested In telecon of 

May 3, 1968 between H.F. Potthast and K.D. Albert. 


Revised Cycle Performance Program (from Reference one) including: 

I o) Computer tape 

b) Program listing 

c) Sample calculation 

d) Input format/ instruct ions 

e) Sample input cards 

The revised cycle performance from the referenced report Incorporates 
the following modifications: 


I 


The program wo s modified in order to more accurately account 
for the pressure distribution along the multi-staged combustor 
walls. This can nesd be. aeccMwp) I shed by stocking: tho input, 
data for each stag;*, spec i fy irsg thie- combustor area raAia and 
equivalent pressure, and- weighing constant it for each sfen^e-. 
The program vflll execute tho multi-stage problem os a single- 
stage problem by using the results from previous stages. • 



NASA Langley Research Center 
Attention: Mr* K. D. Albert 


SPJD-G 267-OS 10 
Page -2- 


2* The momentum-continuity equation In this program was 
modified to a t h i rd order algebraic equation in order 
to more accurately solve for tho velocity at the com- 
bustor exit. 

3* Because the dynamic pressure distribution Is not linear 
along the multi-stage combustor, the -viscous drag cal- 
culation was modified by using the averaged value of q 
a i£_ each. stage rather than using the averaged value of 
q for the overall combustor os done previously. 

Very truly yours, 

AIRES EARCH MANUFACTURING COMPANY 

A Division of The Garrett Corporation 


R. J. Dunn 

Program Administrator 

Advanced Propulsion Engines Sales 

R,iD/lk 

Enclosure 


CC: D. Argenbrlght - DCA50 

M. Craig - Garrett Sales Hampton 



INPUT DESCRIPTION 


Card 


l 

(914) 


2 

(8FI0.0) 


Field 

1 

2 

3 

4 

5 

6 

7 

8 
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Dcscr I pt i on 

Reference case number "I" for first stage 
M I5" for subsequent stages 

”0" means subsonic combustion; " l" means 
supersonic combustion 

Blank 

Blank 

Number of nozzle stations for which data 
are printed out 

Number (0 to 25) of nozzle pressure ratios 
for schedule of assigned nozzle pressure 
ratios 

Blank 

"0" or "1" means normal printout 

"2 M chemical equilibrium calculation 

"3" thermodynamic data for species 

"4" iteration for combustor exit velocity 
and skin f r iction coef f icient 

"5" means both ,f 3" and ''4" 

"6" means all "2", "3" and "4" 

Number of cases (l to 15) 


| "0" means specified combustor pressure ratio 

*'-l" means specified combustor exit Mach No. 
"I" means specified combustor area ratio 

2 A e /A c , ratio of nozzle exit area to cowl 

area. Blank if noizle pressure ratio is 
spcci f i ed 

3 A 2 /A c , ratio of combustor inlet area to 
cowl area. Blank if combustor inlet Mach 
number is specified 



INPUT 


3 


JESCRIPTION (Continued) 


Description 

combustor inlet Mach number. Blank 

if A2 /A c is specified 
« 

P e /P 0 , ratio of nozzle exit pressure to 

ambient pressure. Blank if hi/A c is 
specified t_.» • A* A w 

■ ( 

Cos 0 j 0 is the fuel injection angle with 
respect to the horizontal axis. 

(e.g., "0.” for normal injection) 

A q , cowl area, sq ft 

Av V /A 2, ratio of combustor wetted area to 
combustor inlet area 

Kj, constant for evaluating the average 
pressure acting on the first stage 
combustor walls, p = (P2 + Kj P3 )/( I + Kj ) 

Ambient pressure, atm 

Ambient temperature, °K 

Flight Mach number 

Blank 

Additive drag coefficient 

\v in ' e.t kinetic energy efficiency, or 

(by pre-fixing a minus sign) the enthalpy 
rise between freestream pressure, entropy 
conditions and conditions corresponding to 
freestream pressure and combustor- i n let 
entropy 

A 0 /A c , ratio of inlet capture area to 

cowl area (i.e., mass flow ratio) 

4- — ?. tic/A?,efp #b-f; Wmi-6 

A e , geom/Ac, eff., ratio of geometric nozzle 
exit area to effective nozzle ex ; t area (to 
account for boundary layer growth) 



INPUT DESCRIPTION (Continued) 


i Card 

U. 

(3F10.0) 


7 

(8FIO.O) 


8 

(8F10.0) 


Field 

4 


5 


1 

2 


3 

4 

5 

6 


i 



Descript i on 

C s , nozzle velocity coefficient, ratio 
of I V3C the Ivac calculated with no 
nozzle losses 

Temperature of the fuel at fuel-injection 
nozzle, °K 

<£, equivalence ratio 

77c, combustor chemical efficiency, ratio 

(by weight) of reacted fuel to fuel 
suppl i ed 


Vf, fuel inlet velocity, ft per sec. 


Temperature of fuel corresponding to 
storage conditions . f '/\ 

r' )’ 

Assigned combustor area ratio, /pressure 
( rat ip,; or exit Mach number, as applicable 



i . i 


Combustor skin-friction coefficient 
(a guessed value for multi-stage run) 

l< 2 , constant determining the effective 

pressure acting on the second stage 
combustor wal 1 


The ratio of the first stage exit area 
to the first stage combustor entrance area 



The ratio of the first stage exit pressure 
to the first stage combustor entrance pressure 


Skin-f'r ict ion drag force acting on the 
first stage, Ib^ 

Ratio of the first stage combustor wetted 
area to the first stage combustor entrance 
area 


(Blank for third stage run) 

Second stage combustor skin-friction 
coefficient 

(Blank for third stage run) 
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INPUT DESCRIPTION (Continued) 


Card Field Descrlpt ion 

7 The dynamic head of the flow at the exit 

of the first stage, psia 

(Blank for third stage run) 

9 1 
(8FI0.0) 

2 


The ratio of the second stage exit pressure 
to the first stage combustor entrance pressure 

Skin-friction drag force acting on the 
second stage combustor, lbf 

5 Ratio of the first and second stage 
combustor wetted area to the first stage 

■ combustor entrance area 

6 Third stage combustor skin-friction 
coefficient 

7 The dynamic head of the flow at the exit 
of the second stage, psia 


3 

4 


K 3 , constant determining the effective 

pressure acting on the third stage 
combustor wall 

The ratio of the second stage exit area 
to the first stage combustor entrance area 


^4 

At 


Note: 

Card 8 and 9 are blank for first stage run 
Card 9 is blank for second stage run 



